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CHAPTER I 


INTRODUCTION 

Composites are becoming more widely used as primary structural 
materials, hence the concern for predicting reliability is- also increasing 
As fiber-reinforced composite materials make the transition from labora- 
tory test coupons and non-critical components to major production items, 
methods for predicting the fracture and strength characteristics are 
necessarily drawing more attention. The question of attempting to under- 
stand the nature of fracture and damage growth in laminated composites 
can be approached in several different ways. One may analyze the laminate 
as being composed of a homogeneous orthotropic material (or of homogeneous 
orthotropic layers) and use numerical methods such as finite elements 
to obtain results for the response of particular structures or even to 
investigate local effects near geometric or load discontinuities. The 
use of an equivalent homogeneous material and of numerical methods 
to investigate local effects is not always satisfactory [1]. It is often 
the multiphase nature of a composite that gives it its unique properties, 
especially near the damaged regions. To replace the actual material 
by an equivalent homogeneous material may well remove one of the most 
important characteristics of the material. 

Another approach is to use empirical models (descriptive models) 
which generally do not account for the physics of the particular materials 
but contain parameters that are adjusted to gain results agreeing with 
experiments. Such descriptive models have little use in understanding 
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damage growth and local effects, but can be valuable in developing 
design guidelines for structural components. 

By contrast, the present study is concerned with developing rela- 
tively simple mathematical models that contain the important physical 
and geometric properties of the composite (predictive models) such that 
it is possible to investigate the fundamental behavior of a laminate 
in terms of the various material properties of each lamina. One of 
the main difficulties with such modeling is selecting significant 
properties to have a model complete enough to be reasonably accurate. 

At the same time it should be simple to solve. Zweben gives an excel- 
lent discussion of this "materials modeling approach" in [2]. 

Some of the first work in modeling a unidirectional composite in 
this manner was done by Hedgepeth [3], where no damage other than an 
initial transverse notch was considered. The study was extended by 
Hedgepeth and Van Dyke [4] for the special case of one broken fiber 
with longitudinal splitting. The extension to more than one broken fiber 
with longitudinal damage could not be developed conveniently by 
influence functions as in [3] and [4] because the broken fiber adjacent 
to the damage region was not typical of any of the remaining fibers. 

Goree and Gross [5] used Fourier transforms to modify the solution to 
account for an arbitrary number of broken fibers as well as for longi- 
tudinal matrix damage to include both yielding and splitting initiating 
at the notch tip between the last broken fiber and the first unbroken 
fiber. Goree, Dharani and Jones [6] added constraint layers to the main 
lamina to account for either a misalignment of fibers in a multi-ply 
unidirectional laminate, or for the presence of angle plies which give 



3 


support to the unidirectional ply. The constraint layer was taken as 
being fully bonded to the unidirectional lamina at all times in the 
above model. 

The analysis presented here is an attempt to extend the above dis- 
cussed model to include the effect of debonding between the notched 
unidirectional lamina and the surface constraint plies. The laminate 
is modeled as a two dimensional region of a unidirectional lamina 
with symmetrically located surface constraint layers whose fibers make 
an angle theta with the unidirectional ply, (Figure 1). 

In the vicinity of a notch in a laminate, the broken fibers exert 
longitudinal shear stresses in the matrix which are transferred to the 
nearest unbroken fibers. The shear-lag assumption [2, 7, 8] is used 
for this shear transfer between fibers in the unidirectional lamina. 

An equivalent expression is used for the corresponding effect of the 
constraint layers. Since the shear transfer does not depend on the 
transverse displacements, it uncouples the longitudinal and transverse 
equilibrium equations. Hence, the longitudinal displacements can be 
calculated without solving the transverse displacement equations. 

The fibers support all of the longitudinal stress in the unidirec- 
tional lamina because the longitudinal extensional modulus of the 
fibers is assumed to be much larger than that of the matrix. Debonding 
between plies in the vicinity of the crack is considered to be of 
finite width and extends to infinity in the longitudinal direction. 

The extreme fibers of the debond zone are assumed to be attached to each 
other across the debonded region by a spring of stiffness depending on 
the constraint layer properties and on the width of the debonded zone. 
That is, the layer debonds but is still connected to either side of 
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Figure 1. Unidirectional Lamina with Broken Fibers, Surface 
Constraint Layers and Debonding. 
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the region and carries load due to the longitudinal displacement of 
the extreme fibers of the debonded zone. This can be more clearly seen 
by referring to Figure 1. 

As an initial investigation, the basic mechanism of crack growth 
is limited to a model containing broken fibers only. Subsequently, a 
model is next developed to account for additional longitudinal damage 
parallel to fibers in the monolayer, (Figure 2). Splitting and yield- 
ing of the matrix is assumed to initiate at the notch tip and to pro- 
gress longitudinally between the last broken fiber and the first 
unbroken fiber. 

There are three different zones in the model. 

(i) unidirectional ply with bonded constraint layers, fiber 
numbers (0) to (N-l) and (M+l) to (°°). 

(ii) unidirectional ply with debonded constraint layers, fiber 
numbers (N+l) to (M-l), and 

(iii) intermediate fibers (N) and (M) of the following ply. 

An equilibrium equation is written for each fiber using the basic 
stress-displacement relations given by Hooke's law and the shear-lag 
assumption. A description of the solution to these equations is given 
in the following sections. The stresses and displacements are determined 
as a function of number of broken fibers, constraint layer parameters 
and debonded zone size. The results are compared to the corresponding 


fully bonded cases. 
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Figure 2. Unidirectional Lamina with Broken Fibers, Longitudinal 
Matrix Damage, Surface Constraint Layers and Debonding. 


longitudinal damage 
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CHAPTER II 
FORMULATION 

Two Dimensional Shear-Lag Model with 
Broken Fibers, Surface Constraint 
Layers and Debonding 

A unidirectional array of parallel fibers with surface constraint 
layers, debonding and an arbitrary number of broken fibers is shown in 
Figure 1, Debonding is assumed to exist from the last broken fiber (L) 
to an unbroken fiber (M) of the unidirectional lamina. The constraint 
layers are intended to represent adjacent layers of a unidirectional 
lamina. They are assumed to be placed symmetrically about the uni- 
directional lamina to give a laminate with no bending. The broken 
fibers are assumed to occur along the x-axis and, since the loading 
is symmetric, only the first quadrant of the laminate is considered 
in the analysis. The basis analysis and assumptions are the same as 
in [6], however, in order to clearly indicate the modifications needed 
to account for surface debonding, it is necessary to repeat the basic 
formulation. 

The fibers are taken to be of much greater stiffness and strength 

than the matrix and the longitudinal load is therefore assumed to be 

carried by fibers only. Load is transferred between fibers by shear 

stresses as given by the classical shear-lag assumption. The axial 

fiber stress a (y) and the matrix shear stresses, t (y) and T ! (y), 
n n n 

are then given by the simple relations: 
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dv 

a (y) = E -j— • » 
n w F dy 


M 


T n (y> = T [v n (y)-v n _ 1 (y)]. 


x;(y ) 


G m 

M 


= [v n (y)- v n_i(y)]» 


where 


(1) 

( 2 ) 
( 3 ) 


v^(y) - axial displacement of the fiber (n) at the location (y) , 
Ep = Young’s modulus of the fiber, 
t - thickness of the unidirectional ply, 
t 1 = thickness of the constraint plies. 


The stiffnesses G^J h and G’/h T must account for interaction 
between fibers [5, 7, 8]. G^. and G^ are typically not the shear moduli 
for the "neat" matrix nor are h and h f necessarily fiber center-line 
distances. The ratios G^/h and G^/h 1 are equivalent stiffnesses and 
are assumed to be material constants and depend only on the fiber and 
matrix properties, the fiber volume fraction, orientation of plies, 
and not on the size of the damage region. Only for large spacing can 
and h be expected to approach the "neat' 1 matrix and center-line 
values . 

By the virtue of the shear-lag assumption the longitudinal and 
transverse equilibrium equations become uncoupled and the longitudinal 
displacement and stress in the fibers as well as the matrix shear 
stress can be obtained without solving the transverse equilibrium 
equations. Therefore, only the equilibrium equations in the longi- 
tudinal direction will be considered in the following discussion. 
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As pointed out earlier and in [5], G„/h and G'/h' are to be de- 
termine experimentally. For example in [9] and [10], it is shown that 
the shear stress becomes larger as the fiber spacing decreases, that 
is, (0(l//d)) for rigid fibers where ’d' is the minimum distance be- 
tween the fibers. Local failure may occur at critical points through 
the thickness in advance of laminate splitting which would give an 
apparent shear stiffness considerably different from that of the matrix 
alone. 

The debonded fibers, (N) and (M) are considered to be connected 
by springs due to the presence of the angle-plies of the constraint 
layers. The springs are assumed to have a linear force-displacement 
relation and the stiffness (k) per unit area for a particular laminate 
to decrease, proportional to the length of the spring. 

With reference to the free-body diagrams, Figures 3 through 5, of 
the elements for different ranges of fibers, the equilibrium equations 
are 


d °F 


+ (t| 


- T 


t dy ‘ Vl |n+1 
for n = 0, 1, 2, . . . ,N-1, M+l,..., 


) + (T'| ) -t- = 0 

n [n-KL |n t 


T^r + (t | n+l ~ T 


n 


) = 0 


for n = N+l, N+2, . . . ,M-2, M-l, 


where 


Ap = area of fiber. 


For fiber N 


d °F , | 1 . , kt' . , r t f 

t W + <t |n+i- t |n ) +_ r (v m-V- t |nT = 


0 . 


(4) 


(5) 


( 6 ) 





NOTE: 



Diagram shown as above for clarity only. In the actual laminate, 
the constraint layers are symmetric with a layer of thickness 
t ' / 2 placed on each side of the unidirectional ply. 


Figure 4. Free-body Diagram for a Typical Element of the Fully Bonded 
Zone (Fibers No. (0) to (N-l) and (M+l) to («>)). 





NOTE: Diagram shown as above for clarity only. In the actual laminate, 

the constraint layers are symmetric with a layer of thickness 
t'/2 placed on each side of the unidirectional ply. 


Figure 5. Free-body Diagram for a Typical Element of the Intermediate 
Zone (Fiber No. (N)). 
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For fiber M 


tx % 

t dy 


+ (t 


M+l 



kt T 

t 


(v m" v n ) + t ' 


M+l 



(7) 


Using the stress-displacement relations (1) , (2) and (3) , in the 
above equilibrium equations, the following set of difference-differen- 
tial equations is obtained: 


W d v n 


G t , 2 + (1 + C R ) ( Vr 2 v n + Vl } = ° 

M dy 

for n = 0, 1, 2 , . . . ,N-1 , M+l,..., 


( 8 ) 


A F E F h d V n 

r 7T + (v - - 2v + V ) = 0 

G t ,2 n+1 n n-1' 

M dy 

for n = N+1, N+2 , . . . ,M-2 , M-l. 

For fiber N 


(9) 


A F E F^ d V N 

G M t dy 2 + (v N+l “ 2v N + V N-1^ + C R2 (v M _ V N ) " C R (v N~ V N-1 )=0, (10) 


For fiber M 


A F E F h d V M , 

G M t dy 2 + (V M+1 2V M + V M-1 ) C R2 (v M V +C R (v M+1 _ V M )=0, (11) 


where 

C R2 = Tc^hT "t”’ 

( 12 ) 

(G M /h,) t* 

Cr (G M /h) t • 

The constraint layer provides additional longitudinal stiffness to 

the unidirectional ply, the effect being given by the constant C_. The 

R 

debonding effect is represented by the second constraint ratio C 0 . 

R2 
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To match the differential equations for the bonded case [9], C r 

must reduce to when M-N=l. The only varying parameter for 

in a particular laminate is the width of the debond zone, (M-N). 

Since represents a linear spring and loses its stiffness (k) per 
RZ 

unit area proportional to its length then C R2 “ M-N ' 


Hence, C 


R2 M-N 


Noting the coefficient of the second derivative in Equations ( 8 ) 
through ( 12 ), the following changes of variables as suggested in [ 6 ] 
are made. 


A F E F h 


a = a a = E„ - 5 — * 
n 00 n F dy 


V 1 

V = Mw v 


where 


a = applied remote stress, 

CO 

and n, a and V (n) are non-dimensional. 

’ n n 

Algebraic manipulation of Equations (3) and (14) gives 


a = a - 3 — > 

n 00 dy 


g m a f . . 

(v n- v n-l>- 

r 


The resulting equations in non-dimensional form are 


T + Ci + c R ) -o 


for n — 0, 1, 2 , . . . ,N 1, M - ! - 1 , . . • , 
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d 2 V 


+ (V , , - 2V + V , ) = 0 
2 v n+1 n n-1' 


dn 


for n = N+1, N+2,...,M-2, M-l. 


For fiber N 


d 2 V. 


dn 


For Fiber M 


r + <Vi' 2 V v n-i )+c r2 ( VV- c r ( V v n-i> ' °- 


d 2 V. 


TT + ( Vi- 2V M + W- c R2 (v M-V- c R<Vi-V ■ °- 

dn 

These equations can be written as follows, where the left-hand 
side is the same in each equation* 


d 2 v 


dn 


f + < 1 +C R><Vl- 2V n + Vl> ‘ 0 


for n = 0, 1, 2,.,.,N-1, M+l, . . . , 

T* ( ltC R>‘Vl-Wl) - V V n+l- 2V n + V n-l> 


d 2 V 


dn 


for n = N+1, N+2 , . . . ,M-2 , M-l. 


For fiber N 


d 2 V, 


(18) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


dn 


For fiber M 


2 + ^ 1 + C R^ V N+1 2V N + V N-1^ Wl V- C R2< V M-V- < ‘ 23 ^ 


d 2 V, 


M 

— f + (1 + c r )(v m+i- 2v m +v m-i> =- c R (v M- v M-i ) + c R2 (v M-V- (24) 
dn 

These difference-differential equations may be reduced to differen- 
tial equations as in [6] by introducing a new function V(n»9) defined 
as 
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v (n) 

V(n,0) = — 2 + 2 V n^ n ^ cos 

n=l 


from which 


V (n) = — / V(n>9) cos (n0)d0. (26) 

n it Q 

Making use of orthogonality of the circular functions. Equations 
(21) through (24) are then written as one equation, valid for all q and 


— / f - 2(1+ C R ) [1 - cos(0) ] v'l cos (n0)d0 

77 o \ dq K J 


7T r h 
0 U=N 


G^(n) cos (£0) cos (n9) > d0, 


V’ “ - C R (V t+l- 2V t + Vl> 


for f. = N+l, . . . ,M-1, 


G n (0) - -V V N + 1-V +C R2 <V M-V> 


C R< V M- V H-X>- C R2 <V M-V’ 


G 0 (n) - - y IWV- c r2< v m- v o»- 


The equation is of the form 


— / F(n,0) cos (n9)d0 = 0 for all n and n , ( 

17 0 

and as F(q,0) is even valued in 0, if the integral is to vanish for 
all n> the function F(q,0) then must be zero. The single equation 


specifying V(q,0) is then 

2- M 

- 6 2 V = - z g (n) Cos Ue) 
dn £=N 


where 


= 2(1 + C R ) [1- cos (0)]. 


( 33 ) 
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The solution to the problem of vanishing stresses and displace- 
ments at infinity and uniform axial compression on the crack surface 
will now be sought. The complete solution will be obtained by adding 
the results corresponding to uniform axial stress and no broken 
fibers to the solution of the following problem. The appropriate 
boundary conditions are: 

dv n (n) 

V n (n) = 0 , — ^ — =0 as n + °° for all fibers, (34) 

dv n (T l) 

— dn — = “ ~ 1 at n = 0 for all broken fibers, (35) 


V (ti) =0 at n = 0 for all unbroken fibers. (36) 

The complete solution to Equation (24) , satisfying vanishing stresses 


and displacement at infinity is given by 


M 


V(n,6) = A(0)e (6n) + / j sinh(6(n-t)) 2 G. (t)cos(£0)dt , (37) 

n £=N 

where the function A(9) is yet unknown. The remaining two boundary 
conditions give 
dV (0) 


7 (0) TT r oo/' M \ -i 

~rr — = ~ f -6A(0)+ / <cosh(6t) Z G 0 (t) cos (£0)dt ) cos (n0) I 
n 0 L 0 v £=N * J J 


d0 


= - 1 

for all broken fibers, and 


(38) 


V (0) = \ f Fa(0) - t / /sinh(6t) Z G 0 (t) cos (£0) dtl cos (n0)l 

0 *- ° 0 l £=N* ; J 


d0 


= 0 

for all unbroken fibers . 

Equation (37) can be solved exactly by taking 

i 00 M L 

A(0) — ~r f sinh(6t) Z G„ (t)cos(£0)dt = Z B cos(m0). 
5 0 Jl-N * m=0 m 


(39) 


( 40 ) 
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Eliminating A(0) from Equations (38) and (40), the stress boundary 
condition (35) reduces to 

Tf ✓* oo \ 

— f (-6 E B cos(m9 )+ / e Z G„ (t) cos (£0) dt / cos(n0)d0 (41) 

* 0 l m=0 “ 0 £=N J 


= - 1. 

From Equations (37) and (40) , A(0) can be eliminated to obtain 

V(n,e) in terms of B and G (t) . Recalling the relation between 

m x 

V(n,0) and ^^(n) an expression can be obtained for longitudinal fiber 
displacements as 

V 


= — j (e Z B cos(mO) 


■ i ' {• 


m=0 


m 


+ 


E G 0 (t)cos(£0)dt > l cos(n0)d0 (42) 
2 0 5 £=N 11 J 


where 


^ “ 6 h-t -«(n+t) 

D(6,q,t) = e - e 


(43) 


The longitudinal fiber stress is obtained by differentiating 
Equation (42) with respect to n and is 
dV (n) 


5 n (n) - 


n 


dp 


+ f 
0 


2 77 r L 

- f 1-6 £ 

* 0 V m=0 

r,*. 


B cos(m0)e 
m 


»(6n) 


-6(n+t) 


- p 


where 


E G p ( t) cos (£0)) cos(n0)d0, 
H=N 

f 44^ 


») 


p = 1 for t < T] ? 
p = - 1 for t > n # 

Equations (28) , (29) , (30) , (41) and (42) can be solved for the unknowns 
B m and G £ (t) * 
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Limit Case of an Infinitely Wide 
Debonded Zo ne 

This is an extension of the model developed previously, such that the 
debonded zone is now assumed to extend to infinity. Since by Equation (12) 


r _ k t ' 
R2 (G M /h) t ’ 


ko: M=N’ then C R2 ^ 0 
for a debond zone of infinite width. 

Physically, the spring between fibers at the extremities of the de- 
bonded zone has no stiffness as it has an infinite length. Also, for 
fibers (M) and (M-l) far from the crack tip V M = V . Then using 

C R2*° V M- V M-1 

Equations (21), (22) , (23) , (24) reduce to 


U + C R ) (V n+1 - 2V n + V n _ x ) - 0 
for n - 0, 1, 2 , . . . ,N-1, M+l and 


T + ( V n + l- 2V n + V n-l> = 0 


for n = N+l, N+2 , . . . ,M-2 , M-l 
For fiber N 


T- <Vi- 2v H + v n-i> - - C R<Vl-V- 


Following the same technique as before, the single differential 
equation to be solved is then 

d 2 V 2 - N 

— o " 5 V = - £ H 0 (n)cos(£9) 

dq £=0 * 
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where 


V") ' VV 1 - 2v n + Vi ) - 

For fiber (0) 

H 0 (n) ■ C R (V 0- V- 

The boundary conditions (34) , (35) , (36) yield 

if * ^ 00 N \ 

- f (-5 l B cos (m9) + / e" (6t) Z H. (t)cos (£0)dt > cos(n0)d0 

77 0 l m=0 m 0 S.=0 ) 


(49) 


(50) 


(51) 


= - 1 

for £ = 1, 2, 3 , . . . ,N . 

Equations (49) , (50) , (51) , (52) can be solved for the unknowns B m 

and H (t) . The longitudinal displacement is given by 
£ 




e Z B cos(m0) 

* m 
m=0 


7T N > 

1 f Midbil 2 h o (t)cos(£0)dt ) cos(n0)d0, (52) 

2 0 5 1=0 * J 


while the longitudinal fiber stress is given by 


a n (n) 


2 77 f L 

= - f j-6 Z 

77 0 l m=0 


B e cos (m0) 

m 


7T /-6(n+t) 

+ ; o h — p - 


-6 | n-t |\ N 


Z H 0 (t)cos(2,e) > cos(n0)de, 

£=0 J 

( 53 ) 


where 

p = 1 for t < n > 
p = - 1 for t > n • 
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Two Dimensional Shear— Lag Model with 
Broken Fibers, Longitudinal Matrix 
Splitting and Yielding, Surface 
Constraint Layers and Debonding 

The solution developed in Section I will now be extended to in- 
clude longitudinal splitting and yielding of the matrix as shown in 
Figure 2. All the previous assumptions are assumed valid and it is 
only necessary to account for additional damage. It is assumed that 
splitting and yielding of the matrix initiates at the notch tip and 
progresses longitudinally between the last broken fiber and the first 
unbroken fiber. The matrix is assumed to be elastic-perfectly plastic. 
The last broken fiber is considered to be the first debonded fiber. 

All the equations remain the same as (4) , (5) , (6) and (7) for all 
fibers except for fibers (L) and (L+l). The equilibrium equation for 
fiber (N) is 


^F dg F N 
t dy 


y - lr 



kt 1 
t 


(v — V ) — T T 

V M N ; T 


N 



when y £ , 

for fiber (N +1) 


(54) 


^F ^F N+l , 
t dy + T |n+2 + T 0 


< y- Z 2 


> = 0 


(55) 


when y < £ . 

Tq = matrix yield stress, 

^ = length of longitudinal matrix damage at the crack tip, 

= length of longitudinal matrix split at the crack tip. 

The above equations on introduction of the stress-displacement rela- 
tions (1), (2), (3), become, 
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for fiber (N) , when y _< 


A F E F d ^ _ Xq < y -£ 2 > - -£ (^-V^) + < V M -V 


dy 


and for fiber (N+l), when y 


G* ^ T 

M , x t 

- T <v tT v n-i ) T 


- 0, 


(56) 


A F E F d N+l , G M 

_ IT vv N+2 _ V N+1' ' l 0 
dy 


+ "iT K-u-o - v M4-i ) + T n < y - £ 2 > = °* 


(57) 


Equations (56) and (57) are re-written as 
iLE h d 2 v 

-H- — r + (1 + V <v n+i - 2v n + V i> 

H dy 

= + (1+ C R ) (v N+1 - v„)-C r 2 (v m - v H ) + T 0 < y - t 2 > , (58) 


G M t 


d v, 


N+l 


dy 


+ (1 + C R ) (v n+2 2v n +i v N > 

- (1 + C R ) (v N+1 - v N ) + C r (v r+2 


V N+1> " T 0 


<y- 


%2> . (59) 


l 


1 



a , 


a 


2 





T 


0 



( 60 ) 


The resulting overall non-dimensional equilibrium equations are 
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d 2 V. 


dri 


# + < 1 + V <V n+l- 2V n + V n-l> ‘ 0 


(61) 


for n = 0, 1, 2,...,N-1, M+l,.... 
For fiber (N) , when n <_ a 


d 2 V 


N 


dn 


+ (1 + c r )(v n +1 - 2v h +v n-i> 


- < 1 + C R> <V N+l-V- C R2< V M -V +? 0 <n - 6> - 


(62) 


For fiber (N) , when n > a 


d 2 V. 


N 


2 + (1 + C R )(V N+1 2V N + V N-1 ) C R (V N+1 V C R2 (V M _ V' (63) 


dn 


For fiber (N+l) , when n < a 


d 2 V, 


N+l 


2 + C 1 + C R ) ( v n+2 2V N+1 V N^ 

dn 

= - d + c R ) (V N+1 - v N ) + c r (v n+2 - v N+1 ) - x 0 < n - B >. (64) 

For n = N+2, N+3,..*,M-1, and N+l, when rj > a 


d 2 V 


dn 


f* < 1 + C R )( Vl- 2V n + Vl> -Vrt’W 


(65) 


The differential equation to be solved is the same as Equation (33) 
while G £ (n) is given by 

G N (n) = -( 1+c R )( v N+ i _ V +C R2 (V M _ V “*0 <n_ e> f0r n - a ’ 


C R ( ‘ V N+1 V N^ + C R2^ V M _V N^ 


for n > a , (66) 


°N+1 ^ ) ( ' 1+C R^ V N+1 V N^ C R^ V N+2 V N+P + T 0 < n g> 


for n , 


S' 1 ') ■ C R<Vl- 2 V t+Vl> 

for £ = N+2,...,M-1 and N+l for n > a. 


(67) 

( 68 ) 
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Since G^(q) should match at n = a for £ = N+l and N, 


V N (°0 - v N+1 ( a ) 



(69) 


Equations (41) , (42) , (66) , (67) , (68) and (69) can be solved for the 
unknowns and G^(t) . The expressions for displacements and stresses 
remain as given by Equations (42) and (44) . 
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CHAPTER III 
SOLUTION TECHNIQUE 


In all the problems discussed so far, the solution reduces to a 
series equation coupled with one or more Fredholm integral equations 
of the second kind. Since there is no exact form of solution avail- 
able to solve such a system of equations, a computer program is de- 
veloped by modifying a numerical procedure given in [12]. The tech- 
nique makes use of a method by Riez [13] to solve a linear Fredholm 
integral equation of the second kind defined within a semi- infinite 
interval of integration. A given integral over a semi- infinite in- 
terval may be approximated by the Gauss-Laguerre quadrature rule as 
co k -x i 

f f (x) dx = E w.e f (x. ), (70) 

0 i=l 1 

where x. is the i 1 -^ 1 zero of the Laguerre polynomial, L,(x), and w. 

rC X 

is the corresponding weight function given by 


x. 

X 

w = 2 • 

[(k+l)L k+1 (x.)] Z 


The Laguerre polynomial L^(x) is given by 


(71) 


L fc (x) ■ 


x ,k, k -x x 
e d (x e ) 

A k 

dx 


(72) 


Since the form of the equations for each solution is the same for 
the three particular cases discussed above, the development and appli- 
cation of the numerical procedure can be demonstrated without any loss 
of generality by taking equations corresponding to one of the solutions, 
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Consider the problem of an infinite debonded zone with a con- 
straint layer given by the equations 


V") - Wi- 2 v t + Vi> 


for Z - 1, 2,...,N-1, N, 


V n) = C R (V 1- V’ 


2 "f L 
£ / | -6 E 

17 0 V m=0 


B cos(m0) 
m 


N 


+ / e E H^(t) cos (£0) dt } cos (n6)d0 = 


0 £=0 
for n = 0, 1,...,L-1, L, 


} 


- 1 , 


7 (n) = f /(. 

ol 


e E B cos(m0) 

r , m 

m=0 


+ i e H. (t) cos(J-0)dt'l cos(n0)d0. 

2 o 6 a=o 11 ; 


Substituting the expression for V n (ri) in (49) and (50) yields 


H„(n) = c t 


[?; 


e Z B cos(m0) 

m=0 


X {cos(U+l)0) - 2 cos (£0) + cos ((A- l)0)}d0 


N 


H 0 ( n ) = C R 


1 x " r D(t,n 1 . 0 £ H (t)cos(ne) 

" 0 0 5 n=0 n 

X {cos ( (£ + 1) 0) - 2 cos (£0) + cos ( (£ - 1) 0) }d0dt 
— / e E B cos(m0) {cos(0) - l}d0 

Tf ^ r\ ^ 


] 


0 

77 00 


m=0 


N 


+ -/ f - D -- <S - ? - ^ - > — E H (t)cos(ne) {cos(0) - l}dtd0 
17 0 0 n=0 n 


( 49 ) 


(50) 


(51) 


(52) 


,(73) 


( 74 ) 
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The integral over the spatial variable ’t' in the Equations (51), (73) 

and (74) can now be replaced by the series in Equation (70) to yield 
L 2C nr -(<5n ) ( 

E / e cos(m0) < cos ((2 + 1)0) - 2 cos (20) 

m=0 W 0 L l 


+ cos ( (2 - 1)0) 


}•-] 

E f Kj _ 1 (n. ,t.) - 2Kp(q. ,t.) 

[ it j =1 ^ ^ i J 2 V ‘l’ 3 

+ w t ^ apt.) 

- [-p r x f z^ +1 fn i ,c,) - 2K*|(n 1 ,t ) 

L * „»0 j-1 1 n 1 J " 1 J 

" n t. 

+ K*- 1 (n i ,t.)J e J 


d0 


w H (t.) 
t. n y 
2 


= o, (75) 


for 2 = 1,2, .. . ,N-1, N, 


L 2C it 

e ~r f 

m=0 0 


• -(«n,) 1 

e cos (m0) {cos (0) - l}B m d0 


[hj-^r .f l ( K o (n i’ t j ) ■ K o u i> t j ) } ^ \ H o (t j : 

rs n k r i o fc i 

— EE <K X (n.,t.) - K U (n.,t.) / e ^ w H(t.] 
L* n=l j-1 l n 1 J n 1 J J n J 


= 0, (76) 


and 


L 9 7T 

E — / [- 6 cos (m0) cos (n0)B ] d0 
tt in 

m=0 0 


N k „ TTr-(5t.) t. 

E E — / (e J cos (20) cos (n9) \ e J w G^(t.)d0=-1 

=0 j=l U 0 f J j 


+ 


(77) 


for n 1 » L j 

where 


K m (n.,t.) = / D(6,n. ,t )cos(m0)cos(n0)d0. 
n 1 J Q 1 J 


( 78 ) 
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The solution therefore reduces to solving a system of linear 

algebraic equations for the constants B and for explicit values of 

m 

the functions H^(t) at specified quadrature points. The above system 
of equations is solved by the method of Gauss-elimination with partial 
pivoting. 

The above procedure is followed in all the three cases; however, 
for the instance of longitudinal matrix damage, an additional step is 
required. The presence of longitudinal damage has a very significant 
effect on displacement in the damage region and hence on G^(t). If 
the semi-infinite integrals are expressed as a series expansion given 
by (70) , the number of quadrature points lying in the damage region is 
not sufficient to represent the displacement function accurately in 
that region. For example, out of 35 quadrature points, only five lie 
between zero and two, a typical value for alpha (a). So in order to 
approximate the integral more accurately, more points are required in 
the region from zero to alpha (a) . The following procedure is used. 

The integral is split as sum of two integrals 

00 a oo a 00 

f f(x)d x = / f (x) dx + f f(x)dx = / f(x)dx + / f(x+a)dx. (79) 

0 0 a o 0 


The finite interval integrals can be approximated by Gaussian integra- 
tion [14] as 


f f(x)dx = -| £ w* f(y*), 

0 j=l 3 3 


where 

y J “ (f + i )x*. 

x* is the j th zero of the Legendre polynomial P n (x) , 
and w* is the corresponding weight function given by 


(80) 


( 81 ) 
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2 

d-xf) 


[p;(x*)] 2 . 


The Legendre polynomial is: 


iH 


P (x) = 
n 


O 11 t j n 

2 n! dx 


(x 2 -l) n . 


( 82 ) 


( 83 ) 


Physically it would be more direct to specify the applied stress, 
number of broken fibers, L, and determine the damage zone a and g de- 
pending on given yielding and splitting conditions . As a and g appear 
in the limits of the integrals, this is not convenient computationally 
and it is easier to specify the number of broken fibers, N, and the 
damage zone parameters a and g and to compute the required applied 
stress a . 
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CHAPTER IV 
RESULTS 

Debonding with no Longitudinal Matrix 
Damage in the Unidirectional Ply 

The effect of the width of the debonded zone was of particular 
significance in this study. Results are given in Figure 6 for various 
numbers of broken fibers with a constant constraint ratio. Debonding 
was assumed to start at the last broken fiber and extend longitudinally 
to infinity , (Figure 2). The critical fiber is defined as the fiber 
which has the maximum stress. The stress in the critical crack-tip fiber 
decreased initially for a small debonded zone, but subsequently increased 
with an increase in the width of the debonded zone. In fact, the stress 
concentration in the limit case of an infinitely wide debonded zone, 
was more than that of the bonded case, (Table I). 

One result of particular significance is that the maximum decrease 
in the stress in the critical fiber occurs for a small debonded zone 
and is essentially independent of the initial crack length. Figure 6 
shows that a debonded zone width of two fibers spacings results in the 
largest decrease in the maximum fiber stress for five, seven, nine and 
twenty-one broken fibers. 

Debonding acts like a constraint (C R 2) between the last broken fiber 
(N) and last debonded fiber (M) , resulting in a redistribution of stresses 
in the vicinity of the crack tip, hence decreasing the stress in the 
critical fiber. A higher constraint ratio resulted in a larger drop 
in critical stresses as shown in Figure 7 but gave higher critical 
stresses in the limit case Figure 8. 



^F/N 7 ^QO 


Constraint Ratio = 0.5 

Debonded Zone Starts from 
Last Broken Fiber. 

Crack Length = No. of Broken Fibers 
x Fiber Spacing 



debonded zone width (M-N) , fibers 

Figure 6. Maximum Fiber Stress as a Function of Debonded Zone Width 
and Crack Length. 




debonded zone width (M-N), fibers 


Figure 7. Maximum Fiber Stress as a Function of Debonded Zone Width 
and Constraint Ratio. 
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Figure 8. Maximum Fiber Stress as a Function of Crack Length 
and Constraint Ratio. 
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Table I. Maximum Stress Conconcentration vs. Width of Debond Zone 


Number of Broken Fibers = 7 

First Debonded Fiber (N) = 3 

Constraint Ratio (Cr) = 0.5 


Width of Debonded Zone 

Maximum Stress Concentration 

(M-N) fibers 

K 

1 

2.5461 

2 

2.3258 

3 

2.3392 

4 

2.3601 

CO 

2.5813 

Debonding with 

Longitudinal Matrix 

Damage in the 

Unidirectional Ply 


The effects of debonding accompanied with longitudinal matrix splitting 
and yielding at the crack tip are indicated in Figures 9 and 10 where 
some typical results were obtained for seven broken fibers. A two/one 
split strain to yield strain condition is assumed. This ratio was selected 
for comparison with the results of [5] and is approximately equal to 
that for brittle epoxy. A debonded zone of two fiber widths starting 
at the last broken fiber is assumed. 

The maximum fiber stress, normalized by a laminate constant 


To = x 0 / is plotted against the normalized applied stress. 

a J G m A f 

Figures 9 and 10 give results for Cr = 0.5 and 1.0, respectively. 

The results are plotted for a monolayer having four different combina- 
tions of constraint and/or damage as given below. 



Transverse 

Notch 

Constraint 

Layer 

Matrix 

Damage 

Debonding 

o 

X 




ii) 

X 


X 


iu) 

X 

X 

X 


iv) 

X 

X 

X 

X 





Number of Broken Fibers = 7 


Cr - 0.5< 


yield split ^(without debondin^ 


Hedgepeth's^ 

Solution 



X 


Debonded Zone of Two 
Fibers Starts from 
Last Broken Fiber 


Cr - 0.5 
(with debonding) 


T 0 = T 



Cr — 0.0 



no damage 
(net section) 

Op = CJoo 


® yield 
0 split 


Figure 9. Maximum Fiber Stress for Yielding and Splitting for a 
Constraint Ratio of 0.5 



Figure 10. Maximum Fiber Stress for Yielding and Splitting for a 
Constraint Ratio of 1.0. 
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For case (ii), unidirectional lamina with no constraint, it 
was found in [5] that once the split forms the critically stressed 
fiber unloads and the split length becomes unbounded under a five-to 
ten-percent further increase in applied stress. The fracture reduces 
to an unnotched laminate with the net-section fracture stress being 
independent of the initial crack length. 

For cases (iii) and (iv), the critically stressed fiber does not 
reduce to a net section state but continues to carry load after splitting 
with increasing applied load. However, in the presence of debonding, 
the maximum fiber stress is relieved. 

Hence, in terms of load carrying capacity, the worst case is (i) 
where the monolayer has no damage other than an initial transverse 
notch, while the best case is (ii) where the notched monolayer is not 
constrained but has longitudinal matrix damage. Cases (iii) and (iv) 
lie between the above models where debonding allows the larger load 
carrying capacity. 

For all fully bonded constraint layer cases, the maximum 
fiber stress occurs in the first unbroken fiber at the end of the split 
(y = ) for no (or low) constraint ratio [5] and at the notch 

tip (y - 0) for high constraint ratios [10]. The same behavior occurs 
for debonded cases, but in case of high constraint ratios and high values 
of alpha (a) and beta (6), the maximum fiber stress occurs in the last 
debonded fiber at y = 0. This shows that under the above conditions 
the first unbroken fiber is highly relieved of stresses and can result 
in discontinuous damage of the fibers. The following table quantita- 
tively illustrates the above behavior. 



Table II. 


(1) Bonded Cases : 

i ) C R = 0 . 0 . ; a 

n -* 
rj + 

0.00 

3 


ii) C R = 1.0; a = 
n 

ri ^ 

0.00 


(2) Debonded Cases : 

i ) C R = 1.0; a = 


n j 

0.00 


ii) C R = 1.0; a = 


Location of Maximum Fiber Stress 

Number of Broken Fibers = 7 

Broken Fiber Index (L=N) = 3 (Figure 2) 

Stress Concentration at 

n: = Position: n: = Fiber No. 


= 3.00; 3 = 2.25 



4 

5 


1.527 

1.817 

1.447 

1.441 

1.00; 6 = 0.36 


4 

5 


2.294 

2.008 

mmmm 
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CHAPTER V 

CONCLUSIONS AND RECOMMENDATIONS 


The main aim of this research project was to investigate the effects 
of interlamina debonding between notched unidirectional lamina and con- 
straint plies. The significant results are as follows: 

(i) The maximum stress in the crack tip fiber decreases initially 
for a small debonded zone, but increases as the width of the 
debonded zone is increased. 

(ii) The largest decrease in the maximum fiber stress occurs for 

a debonded zone width of two to three fiber spacings. Further 
this debond zone width is essentially independent of crack 
length. 

(iii) For a small debond zone, the maximum fiber stress in the mono- 
layer decreases as the constraint ratio increases. 

(iv) In the presence of longitudinal matrix crack-tip damage and 
constraint layers, the maximum fiber stress in the unidirec- 
tional ply is further relieved if debonding takes place. 

(v) For a large constraint ratio with longitudinal matrix damage 
at the crack tip, debonding reduces the stress in the first 
unbroken fiber such that the maximum fiber stress can occur 
at the end of the debond zone. 

The following modifications are recommended: 

(i) Extend damage to the whole laminate, that is, account for damage 
in the constraint plies. 

(ii) Represent finite longitudinal debonded zones around the crack. 

(iii) Understand the concept of constraint ratio in terms of geo- 
metric and physical layer properties. Since the constraint 
ratio is not defined in terms of layer properties, it is not 
possible to make quantitative predictions about laminate 
strength. Approximate best fit curves [6] and experimental 
studies [15] can help in determining constraint layer 
parameters to find constraint layer ratios. 
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